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X . Introduction . 

A hydromagnetic shock propagating away from the sun is seen by 
a fixed observer as a discontinuous increase in the plasma density, 
n, proton temperature, T, bulk speed, V, and the magnetic field 
intensity, B. It is possible that a similar shock may propagate 
toward the sun, relative to a frame moving with the solar wind, yet 
move outward relative to a fixed observer, the net outward motion 
being due to convection by the solar wind. In this case, one sees 
the rear of the shock first, so n, T and B are seen to decrease 
discontinuously , while the speed V increases. Such a discontinuity 
is called a reverse shock ,, (e.g 0 , see Razdan et al 0 , 1965). 

This paper presents quantitative evidence for such a shock in the 
solar wind. Using Ness's magnetic field data from three spacecraft, the 
shock is shown to be oriented along the spiral direction. It is also shown 
that the shock separates a high speed wind stream from a region of 
gas that appears to have been compressed and heated by the speed 
stream. 

II. Description of the Discontinuity . 

The discontinuity that we shall identify as a shock is shown 
by the dashed line in Figure 1, which is a plot of plasma and magnetic 
field data from Explorer 34. (For a description of the instruments, 
see Ogilvie et al. , 1968, and Fairfield, 1969) , The signature is 
clearly that of a reverse shock. Table I shows that n, T 
and B decrease by nearly a factor of 2 across the discontinuity. 

Higher time resolution data shows that the decrease in B occurred in 
less than ~5 sec at ~2243 UT on 28 Sept., 1967. The change in the 
plasma parameters seems to occur over several minutes. 
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The gradualness may be illusory because of the relatively 

large fluctuations (see Figure 1). The energy 

9 5 

density, defined by P ■ B /8n + nk(T+T g ) with T e * 1*3x10 °K (see 

Burlaga and Ogilvie, 1970 for a review of measurements of T 0 ) decreases 

by nearly a factor of 2 across the discontinuity, primarily because 

of the change in B. 

A 

The magnetic field direction, B, did not change appreciably 
across the discontinuity. From 1.5 min, averages of components of 

— ♦ A 

B measured before and after the discontinuity, we find B^ * (".49, 

.59, .64), and B 2 ■ (".42, .59, .69) in solar ecliptic coordinates 

A A 

(X toward the sun, Z normal to the ecliptic) • The magnitude of B 
did change across the discontinuity while the direction did not, so 
the field has to be parallel to the surface of the discontinuity 
since [B n l =» 0. 

The flow direction was not measured precisely, but the flow did change 

in less than 3 min from a direction >5.9°west of the sun earth line in 

region 1 to <5.9° W. in region 2 at the time of the discontinuity. 

III. Identification of the Discontinuity as a Reverse Shock. 

Basic Equations . To show that a discontinuity is a reverse 

shock, one must show (1) that it propagates relative to the plasma 

in the sunward direction, (2) that v^ decreases across the shock, and 

(3) that the changes in the fluid parameters satisfy the Rankine- 

Hugoniot equations which are as follows: 

[nv >0 (1) 

n 

[fev n v p - B n Bp)/(4TT)]=0 (2) * 

[p v n 2 + P>0 (3) 

+ Y nkT + B p 2 + B n Bp v p^ 

2 Y*1 P 4np 4Ttpv n 


( 4 ) 



The subscripts n and p denote components of vectors normal and 
parallel to the shock surface, respectively; U is the speed of the 
shock surface (assumed plane) in the direction of its normal, 
relative to a fixed observer; p is the mass density# Our method is 
to use (1) and (2) in the determination of U and the shock normal 
n, and then show that the data in Table I satisfy (3) and (4). 

From (1) and (5), we get 

U " "l V ln ~ n 2 V 2n (6) 

n l ' n 2 

From (1) and (2), v * v„ (assuming B *0), which with (5) gives 

lp2p n 

V, - V. (7) 

Ip 2p 

A 

Let a) be the angle between V ^ and n, and let OH-e be the angle 

A 

between and n. Then 

V lp “ V 1 sin (a * e) ( 8 > 

V = V sin a) (9) 

2p 2 


and (7) gives 


e ~ (^2 ^1 ) tan 

V, 


( 10 ) 


for small e. From (6) one finds that U is given by 


U ~ n,V, (cos a) - e sin cos 


nn 


In the next section it is shown that to determine n we need 
1 1 A 

V s U/cos to , where oj is the angle between n and the radial direction, 
s 

— » i 

If V 2 is nearly radial, as it surely is, we may set u) *U); the consequent 

A 

small error does not contribute significantly to the uncertainty in n. 



Thus, 
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V -s- 
g COS 0) 


n^V^C 1-c tan u>) - n 2 V 2 
n* “ n, 


( 12 ) 


1 2 

The quantities e and u) can be determined by the following method 

of successive approximations. Assume e«o and compute V from (12). 

s 

With this, compute u> from (15), (16), and (17) below, and then get a 
better estimate o^ e from (10). Repeat the process until ® and u> have 
the desired accuracy. 

Surface Orientation . The discontinuity in B was observed minutes 
apart at 3 spacecraft, Explorers 33, 34 and 35, separated by tens of 
earth radii, as shown in Table II. If the shock can be approximated 
by a plane on a scale of the spacecraft separations (<100Rg) , the 
measured time delays can be used to calculate the orientation 
as a function of the radial speed of the surface V g . Let 0 be the 
angle between the normal to the boundary surface and the ecliptic 
plane (0 = o implies a boundary normal to the ecliptic, 0>o, if the 
normal is above the ecliptic). Let cp be the angle between the line 
determined by the intercept of the boundary surface with the ecliptic 
plane aid a line in the ecliptic which is perpendicular to the earth-sun 
line. (cp^45° implies a plane along the spiral direction; 0«o, cp=o 
implies a plane normal to the earth-sun line). Then 

(13) 

(14) 


V e « A + B tan 0 

D 


tan cp ** CV + E tan 0 + D, 
s 

where A a (X^ - - t^Yj) , 

B 5 (ZjY 2 - Z 2 Y i y(t 1 Y 2 - t 2 Y x ), C a -t 2 /Y 2 ,Ds -X 2 /Y 2 and Es Z 2 /Y 2 . 

Here (X,Y,Z) is the position of Explorer 34 (subscript 1) or Explorer 
35 (subscript 2) relative to Explorer 33, in solar ecliptic coordinates, 
and the times t 2 are the time intervals during which the boundary 
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moves from Explorer 33 to Explorer 34 and 35 , respectively. 

The coordinates of Explorer 33 are (23.72, 36*73, 11.78). Using 

the data in Table 2 we derive from (1$ and 04) the following relations 

tan 9 - (2.3±.8)xlO' 3 V # (km/.ec) - (.7±.3) (15) 

tan <p - (2.65±.15)xlO' 3 V.(km/sec) - (.04±.02) (1C) 

8 

The components of the shock normal are 

n ■ -cos u> * -cos 0 cos cp 
x 

n « -cos 0 sin cp (17) 

y Y 

n ■ sin 0 
z 

Results . Using the numbers in Table II and the iteration procedure 

o o 

discussed above, one finds e = 4.8 and u> ■ 45.3 . Equation (11) 

o 

then gives U « 270 km/sec . From (12), (15) and (16) one gets 0 - 10.4 , 
cp * 44.4° which with (17) gives the shock normal n ■ (-.70, -.69, .18). 
Thus, the surface is nearly perpendicular to the ecliptic plane and 
aligned along the spiral direction . The speed of the surface along 
n relative to a fixed frame (U = 270 km/sec) is less than the normal 
component of the solar wind on the sunward side * V 2n ■ V 2 cos u> - 
411 km/sec. This gives an important result: The surface is propagating in 
the sunward direction (along - n) relative to the solar wind, at a 
speed ~ 141 km/sec . The normal flow speed behind the shock is 

v^ n * cos (u>fe)-U » 76 km/sec. Thus, the re l ative normal speed 

diminishes across the discontinuity , as it must if the discontinuity 
is a shock. 

The angle between the normal n and the field B on either side 
of the shock is cos *(n*B) - 87°. Thus, the field is parallel to 
the shock surface , as required by the fact that B^=B 2 . Hence, the 
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assumption B *0, used to derive (7), la valid, Aa a check on the 
n 

validity of (7), note that * V^aln(uH-c) - 414 kra/*?ec and 
V 2 p * V^sin u) ■ 416 km/sec. 

Now let us show that the data satisfy (3). This equation can 
be written 


pv 




(3*) 


2 

where P s nk(T+T.) + B / 8tt, which says that the pressure change across 

© 

the surface must be balanced by a change In the momentum flux normal 

1 l 3 

to the surface. From Table I, p|^ * (8.0:11. 5) xlO ^ergs/cm . Using 

2 2 “10 3 
the values of v^ n and v^ n above, we find pv^ ■ 8,6x10 ergs/cm . 

Thus (3 1 ) is satisfied within the errors . 

Finally, consider (4). Since * v^p and B n ■ 0, this can be 

written 

ln!l? e( -Lss® + V)|! (4*) 

2 1 Y-l P 4np ‘ 2 

The change in v_/2 is 7.1x10 "*( cm/sec) . If T « 1.5xlO J °K as discussed 

U C 

above, and if y * 5/3, we find that the change in the RHS of (4') is 
3 2 

5.5x10 (cm/sec) . Since the uncertainty in the RHS of (4*) is certainly 
>50%, we conclude that (4) is satisfied within the errors . 

We conclude that the discontinuity shown by the dashed line in 


Figure 1 is a reverse shock. 
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V. The Shock at a Boundary between a "Driver-Gag" and a "Driven-Gas" 

Figure 1 shows that the shock precedes a high-speed stream in 

which the flow has a nearly constant value ~585 km/sec, which is 

considerably higher than the most probable speed for the solar wind, 

~385 km/sec. (see Burlaga and Ogilvie, 1970). Figure (1) shows that 

a second high speed stream passed the spacecraft nearly a day later, 

preceded by a bulk speed gradient and hot spot beginning ~18 hours 

after the boundary. Between the second (much weaker) interaction 

region and the shock, the flow is more or less uniform with no 

marked density or bulk speed changes and with the temperature and 

magnetic field intensity varying slowly in a complementary fashion. 

5 Q 

Moreover, the temperature T^*2.6xl0 jCis essentially the same as 
that predicted from by the (Vj T) relation. (Burlaga and Cgilvie, 1970) 
Thus, this flow has the characteristics of an advancing high speed 
stream of fresh, equilibrium, solar wind plasma. 

Ahead of the boundary one sees a rather different kind of flow. 

The bulk speed is steadily increasing from 350 km/sec to 540 km/sec, 
the density and magnetic field intensity are unusually high and show 
several abrupt changes, the temperature is anamalously high (this 
is one of the regions designated as a ’’hot spot" by Burlaga and Ogilvie, 
1970) , and the energy density is increasing rapidly and nonuriiformly. 
Thus, this is characteristic of the interaction regions which commonly 
precede high speed streams (Burlaga et al. 1970, Burlaga and Ogilvie, 
1970) with material presumably undergoing dynamical changes. 

It appears, then, that the boundary represents a transition 
between a uniform "driver" gas (possibly a corotating stream) and the 
dynamical material piled up ahead of this gas. 
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TABLE 1 

PLASMA AMD MAGNETIC FIELD PARAMETERS AT THE BOUNDARY 

V(kro/sec) n(cm T(kllo°K) B(v) PxlO +10 ergs/cm^ 

540±5 10. Oil. 0 450il00 12. Oil. 0 14. Oil. 0 

2 585il0 5. 5i. 5 260il00 8.0i.5 6.0i0.5 
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TABLE II 

RELATIVE SPACECRAFT COORDINATES AND TIMES 



i-1 

1*2 

Xi (r e ) 

8.09 

.24 

y i <V 

-46.57 

-92.13 

z i ( r e> 

-14.91 

-6.04 

(min) 

9.2+, 7 

24.0+1.5 
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Fig. 1 


FIGURE CAPTIONS 

The dashed line indicates the reverse shock discussed 
in the text. The plots show connected traces of plasma 
parameters for spectra measured at 3 minute intervals, 
samples of the magnetic field intensity taken at 3 

2 , 

minute intervals, and the corresponding pressure P=B /8 tt 

5 

+ nk(T+T Q ) with T g * 1.5x10 °K. The width of the lines 
shows the extent of the scatter of the observations. 
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